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Abstract 



The present study involves two semi-infinite homogeneous, isotropic, linearly elastic 
plates perfectly bonded end to end. The plates have dissimilar material properties (e.g. 
copper and aluminum). Incident on the interface axe both the lowest symmetric and the 
lowest antisymmetric Rayleigh-Lamb eigenmodes. The resulting transmitted and reflected 
fields contain contributions from all of the real, imaginary and complex plate modes which 
are active at the given excitation frequency. A study of the distribution of energy among the 
various reflected and transmitted Rayleigh-Lamb eigenmodes over a range of frequencies 
is undertaken for a number of material combinations of the bimaterial plate. 
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I. Introduction 



Much of the motivation behind the development of theoretical response curves for 
structures has to do with the nondestructive characterization of their components, and 
in particular, those undergoing fatigue loading and slow crack growth. Early detection of 
insipient flaws through studies of the scattering of high frequency elastic waves has received 
considerable attention in recent years (Ensminger(1988)). 

The forced response of homogeneous plates to surface and buried sources has been 
reported in numerous papers (Ceranoglu et al (1981), Weaver et al (1983), Va- 
sudevan et al (1985)). Detailed analyses have also been carried out on plane layered 
systems including plates and half-spaces (Kennet (1983), Harkrider (1964), Haskell 
(1953)). These studies have shown that the dynamic response of a structure is primarily 
affected by the nature of the source and by the dispersive wave motion induced by the 
geometry of the structure. Early references regarding infinite and semi-infinite plates can 
be found in the monograph by Miklowitz (1978) and the two volume work of Auld 
(1973). 

Recently, the technique of eigenfunction expansions has been exploited in the papers 
of Gregory and Gladwell (1983,1984). These authors successfully tackled problems re- 
lating to semi-infinite plates by expanding the unknown displacement and stress fields into 
infinite series of eigenmodes appended with unknown coefficients. The coefficients were 
found by employing a biorthogonality relationship (Fraser(1976)) between the eigen- 
modes. 

In this paper, a similar procedure is used to study the elastodynamic response of two 
dissimilar semi-infinite plates welded along their lateral boundaries. Working in the fre- 
quency domain, the spectral response is found by representing the unknown displacements 
and stress fields in each plate as infinite sums of Rayleigh-Lamb modes. Numerical results 
are presented for the case of an incident field composed of both the lowest symmetric 
and the lowest antisymmetric Rayleigh-Lamb modes. The distribution of energy among 
the modes in the reflected and transmitted fields is numerically studied over a range of 
excitation frequencies. 
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II. Formulation and Theory 



Shown in figure (la) are two elastic plates of uniform thickness 2 h welded at x = 0, 
with P-wave speeds 01 , 02 ! S-wave speeds /?i,/?2! densities Pi,p 2 ', and shear moduli pi,P2- 
The problem addressed is essentially a two dimensional, plane strain scattering problem. A 
given incident field propagates towards the interface from left to right, which upon striking 
the interface, generates transmitted and reflected wavefields. The objective of the present 
analysis is to study the modal contributions to these fields. For an infinite homogeneous 
plate the only field present should be the incident field. However, when the two plates 
differ in material properties, waves incident on the interface boundary at x = 0 produce 
reflections which combined with the incident waves yield the total field for x < 0. This 
field must match that of the transmitted field (x > 0) along the interface at x = 0. In all 
of the following equations, the time dependence (e~ ,wt ), has been suppressed. 

For a general type of source the incident field can be expanded in the form, 



(i) 



n = l 



+ E Q-uirksOe'*-' 

n=l 



where U^\y)e lknI is the displacement field of the rt tk symmetric Rayleigh-Lamb mode of 
the plate, and the wave number k n is a root of the symmetric Rayleigh-Lamb frequency 
equation. The boldface lettering of U signifies that it is a two-vector. In a similar fashion, 
the subscript AS refers to the antisymmetric modes of the plate and the wave number q n 
is a root of the antisymmetric Rayleigh-Lamb frequency equation. The expressions for the 
mode shapes as well as that for the frequency equations are stated in Appendix A. It is 
presumed that the homogeneous problem can be solved and thus the coefficients P n and Q n 
are known. The first series in the expression for of equation (1) will excite symmetric 
modes while the second series will excite only antisymmetric modes. For the purposes of 
the present study, all coefficients of the incident field are zero except the first term in each 
series which is identically one. However, the mathematical formulation presented below is 
applicable for an arbitrary source field. 

For x < 0 the reflected field can be expressed as: 



(2) 






+ E s-uSw* 






n=l 



n= 1 



which represents leftward travelling waves (the overbar indicates complex conjugation). 
The coefficients Rn and S n are to be determined from the continuity conditions valid at 
the interface. 

The transmitted field for the region x > 0, takes the form : 



( 3 ) 



U (,l (x,y) = E + E 5nU<;ki/)e'*"', 



n — 1 



n = l 
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where the ’hat 1 symbol is used to indicate that there is a different set of elastic constants 
to be considered for x > 0. 

All of the mode shapes in equations (2-3) satisfy traction free boundary conditions 
on the free surfaces of the plates (y = ±/i). For a perfectly welded interface, the total 
displacement and stress fields in the two plates must be continuous across x = 0. 

The symmetric and antisymmetric portions of the fields are split with continuity being 
imposed on each separately. In what follows, only the determination of the coefficients for 
the symmetric case is given. A similar procedure can be performed to determine coefficients 
for the antisymmetric modes. Enforcing stress and displacement continuity across the 
boundary at x = 0 yields: 



(4a) 



(4b) 



(4c) 



£ 

n = 1 



PnU[ n) + Y, R nU ( z n) = £ RnU^ 



ti = 1 



n= 1 






n= 1 



n= 1 



ti = 1 



mi £ + m. £ RJiV = mi £ rJ'",' 



71 = 1 



71=1 



71=1 



OG CSC oo 

(W) Ml £ C„S<" y ' + Ml Y, R " S -v' = Ml £ RJW- 

n = 1 n=l n=l 

In equations (4a-d) U z and U y are the x and y components of the displacement vector 
Us while S zz and S zy are the corresponding components of the symmetric stress vec- 
tor Ss (given in the appendix). The complex conjugates of the displacement and stress 
components satisfy the following relations: 



(5a) 


r/( n ) _ Tf(n) . 

W X ^ X 7 


u (n) 

y 


= u {n) 

c y 


(5b) 


c(n) _ o( n ). 

°xx ~ °XX 7 


c(^) __ 
^xy 


C(") 

J iy 



Using the above relations, equation (4) is rewritten as: 



(6a) 



£(A, - R *M"‘ = 

71=1 



£ R-V 



(") 



71 = 1 
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(6b) 



oo 



oo 



52(P« + Rn )Ul"> = £ 



OO oo 

(6c) /i, = M 2 5] £»$<? 

n=l n=l 

oo oo 

(6d) h,Y,( P "- R ") S W 

n=l n=l 

The infinite systems of equations (6a-d) must be satisfied for —h < y < h. A procedure 
to eliminate the thickness variable (y) from equations (6a-d) employs a biorthogonality 
relation valid for elastic plates. This relation is stated in Appendix A as equation (9A). 
Equation (6a) is multiplied by y.\ (where the superscript (m) refers to the m th mode) 
and added to equation (6d) multiplied by —Uy m \ The resulting equation is then integrated 
with respect to y from —h to +h. In similar fashion, equation (6b) is multiplied by /Z 2 

equation(6c) is multiplied by —Ui™' and the two are added and integrated over the plate 
thickness. The biorthogonality relation reduces the resulting pair of equations to the form 



( 7 ) 



( Rm ~ Pm)Jm = ^ R n ^r 



n— 1 



( 8 ) 



OO 

^ ] ( P n Rn)L mn = P-mQm 

n-l 



In the above, 



(9a) 



n = f" + ^U^S^dy 

J — h 



(9b) L mn = -K nm . 

In order to cast equations (7) and (8) into tensor form diagonal matrices Ji } and Qij are 
introduced and defined in equations (12a,b). The quantities J m and Q m in equations (7,8) 
are the nonzero diagonal entries defined in equation(12). Notice that whereas J m and Q m 
involve only the properties of the left and right hand side of the plates respectively, the 
coefficient K mn is a coupling term that combines the properties of both materials. 
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Equations (7) and (8) constitute an infinite system of equations for the infinite number 
of unknowns R m and R m . To proceed with the numerical calculations, the infinite series 
is truncated after N terms. Ideally, the value of N should be chosen in such a way that 
the magnitude of the truncated terms is smaller than a specified tolerance. In the present 
case a practical procedure was adopted. If an increase in the value of N was found to to 
have an indiscernible effect on the output, that value of N was chosen as final. In physical 
terms truncation implies that the influence of higher modes is relatively small. Rewriting 
the truncated equations (7-8) in tensor notation yields: 



( 10 ) 



J,jRj — KjjRj ; — JijPj 



(11) L j; R ; — QijRj — — 

where and Q |; are NxN diagonal matrices involving material properties from the left 
and right plates respectively; 



[0 i*j | 

p J\~uys$ + uys£)dy i = jj 

( o * # j 1 

I ^ + W&'Wy < = ; j : 



and K tJ and L t; are NxN full complex matrices given in equations (9a-b). The quantities 
R and R axe N dimensional vectors containing the unknown reflection and transmission 
coefficients respectively. The right hand side of equations (10,11) involve the vector P 
which contains the known incident field coefficients. Letting Z = J U P ; and Z = — L U P ; , 
equations (10) and (11) can be written as a 2Nx2N matrix equation as follows: 



(13) 



J L r 


R 




Z 


L -Q 


R 




Z 



The task of determining the symmetric reflection and transmission coefficients is re- 
duced to solving equation (13), or, as is done in the present work, solving equations (10) 
and (11). A similar process yields the antisymmetric coefficients. 

Having obtained the coefficients, an analysis of the stress and displacement fields in 
either plate can be undertaken. Of particular interest is the fraction of the incident energy 
flux carried away by each real reflected and transmitted mode. Following Gregory and 



Glad well( 1983), the fraction of energy flux carried by each real eigenmode at a particular 
frequency is given by : 



(14a) E, = |/m(J ) )//m(J 1 )III^H 2 , 

for reflected, and 

(14b) £ J = |Jm(Q,)//m(./ 1 )|||jj,|| 2 

for transmitted modes respectively. ’’Edge” modes corresponding to Rayleigh- Lamb eigen- 
values with nonzero imaginary components transmit no energy and therefore their energy 
fractions are set to zero. It must be kept in mind however, that while the imaginary and 
complex modes do not transport energy they affect the magnitude of the coefficients R 
and R and therefore the energy partitioning. 
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III. Numerical Procedure 



In order to set up the matrix equations (10-11) of the previous section it is necessary to 
calculate the roots of the Rayleigh-Lamb frequency equation, from which the mode shapes 
can be determined and ultimately the matrix entries. In the present work a method 
which closely parallels that of Gregory and Glad well (1983) is adopted which tracks 
a prescribed number of roots over a range of frequencies. These root finding routines (for 
both anti-symmetric and symmetric modes) are called intermittently by a main program 
which solves the matrix equations ( 10 - 11 ) and calculates the fractional amount of energy 
transported by each real mode using equations (14a,b). 

Initially the roots of the following two equations: 



(15) sinh(2k n h) ± 2k n h = 0 

are found for a prescribed number of modes. The equation with the ’ plus” sign is associated 
with the symmetric roots of the Rayleigh- Lamb frequency equation as the frequency ap- 
proaches zero, while the equation with the negative sign corresponds to the anti-symmetric 
Rayleigh-Lamb equation. 

As the frequency is increased, updated values of k n axe found by employing a Newton 
iteration scheme when k n is complex and a bisection algorithm if k n is either purely 
imaginary or purely real. When the real or imaginary parts of k n change sign, the value 
of k n becomes purely real or purely imaginary. In such instances, a left hand turn rule has 
been adopted to ensure that modes propagating with the correct group velocity have been 
selected. 

Figures (lb, lc) illustrate the paths taken in tracking a few of the lowest symmetric 
and anti-symmetric modes of a plate of glass which has a Poissons’ ratio of nearly .25. 
Each of the curves have been labelled with a number (signifying the mode) and a letter 
(identifying that portion of the curve as complex, real, or imaginary). The base axis on 
the right is the real axis, while that on the left is the imaginary axis. The vertical axis is 
frequency in megahertz. While the vertical axis is labelled up to 1 megahertz, the curves 
do not have data points beyond .738 megahertz for the antisymmetric roots, and .696 
megahertz for the symmetric roots. In addition, curves extending beyond the boundaries 
of the graph have been truncated. 

Rather than give a detailed explanation and description for each curve, attention will 
be focused on the third symmetric root since it displays as complicated a route as can be 
seen in any of the remaining curves. The third symmetric root (£' 3 ) is initially complex. 
At .275 megahertz k$ intersects the real plane and becomes purely real. At this point, the 
third mode has a negative phase velocity and a positive group velocity. At .291 megahertz, 
which agrees with the theoretical value of the first cutoff, £3 intersects the imaginary plane 
and becomes purely imaginary. It makes a loop in the imaginary plane, intersecting the 
real plane at .336 and becomes purely real (this value corresponds to the second theoretical 
cutoff). 
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Numerical Results and Discussion 



In each numerical simulation left to right travelling symmetric(S) and antisymmet- 
ric AS) fundamental modes with unit energy flux are incident upon the bimaterial plate 
interface. The unknown amplitude coefficients in equations (2-3) are found by solving the 
matrix equations (10-11), from which the energy fractions can be computed using equations 
(14a, b). 

The values of N are fixed at 21 for symmetric and 20 for antisymmetric modes in each 
of the numerical experiments reported below. Greater values of N than those chosen were 
found to produce indiscernible changes in the amounts of energy fractions as displayed in 
the graphs to follow. The choices for N were determined from a sensitivity analysis of the 
output to varying N. Results of this analysis are presented in tables I and II, where the 
two materials are aluminum and steel (corresponding to the first numerical experiment 
described in detail below). The tables include the energy partitions of the antisymmetric 
reflected and transmitted modes at two frequencies (results for the symmetric case show 
equivalent convergence). The higher frequency corresponds to a Rayleigh wavelength of 
about 1/3 of the plate thickness and the lower frequency about twice that. Tabulated are 
the energy fractions for the first four modes found by taking N equal to 2, 6, 10 and 20 in 
the calculations. Notice from the tables that there is not a substantial difference between 
the convergence at 10 and 20 modes with between two to four significant digits remaining 
unchanged. 

It should be noted that in the following table as well as for all of the results presented, 
the sum of the energy fractions equal 1 to at least five places (for symmetric as well as 
antisymmetric cases) for each value of N and at every frequency point for which the energy 
was evaluated. This provides a small measure of confidence in the numerical calculations 
as well as indicating that energy conservation cannot be used in determining convergence 
of the eigenfunction expansion series. 

The first pair of materials considered are Aluminum and Steel with the incident wave 
originating in Aluminum. In all of the examples considered the plate thickness is taken to 
be 1 cm. The elastic properties of Aluminum used are: a P-wave speed of 6.15 Km/sec, an 
S-wave speed of 3.10 Km/sec, and a density of 2.5 gm/cc. Aluminum is soft in comparison 
with Steel for which the P-wave speed, S-wave speed, and density are taken to be: 5.76 
Km/sec, 3.16 Km/sec, and 7.8 gm/cc. In figures 2 and 3 are shown the fraction of incident 
energy carried away by each mode for the symmetric and antisymmetric excitations re- 
spectively. The upper plot in each of the following figures represents the reflected field and 
the lower plot the transmitted field. The upper limit of frequency chosen is high enough 
to include the effect of up to 4 real modes in either plate material. 

A noteworthy feature of the behaviour of the lowest S mode in figure 2 is the sharp 
maximum in the reflected field in Aluminum just beyond the point denoted by Ri and 
the corresponding minimum in transmission just beyond T\. This sharp maximum in 
reflection also marks the beginning of the negative phase velocity branch of the third S 
mode for steel shown by the dashed curve. The next prominent feature for the lowest 
mode occurs near i? 2 ,T 2 which marks the beginning of the cutoff of the third S mode 
for Aluminum. The peak in the dashed curve near i?2 corresponds to the local minimum 
near T 2 in the fundamental mode in transmission. Shortly thereafter the third real mode 
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for Steel becomes activated causing the transmission curve for the first S mode to rise 
again. The next fluctuation in the curves occurs in the region of 7? 3 ,T3. Again, observe 
that the fourth S mode for Aluminum is activated before its counterpart in Steel. Hence 
there is first an increase in the reflection curve for all four modes in Aluminum. Soon 
thereafter the fourth S mode in Steel becomes active causing a rise in transmission with 
a corresponding drop in reflection. This feature is characteristic of ail the curves where a 
maxima in the reflected energy corresponds to a minima in the transmitted curve or vice 
versa since the total energy must add up to the incident flux. A general feature of figure 2 
is that one sees a blip on the curves for either material at the cutoffs of the other material. 
This indicates that as a new energy bearing real mode becomes activated in one of the 
materials there is seen a fluctuation in the energy curve of the other material. This is 
because of a redistribution of energy among the real modes. Similar features are observed 
in the subsequent plots. 

In figure 3 is shown the results for the antisymmetric case. Once again the cutoff for 
the second AS mode in Aluminum occurs before that in Steel. Thus there is first a rise 
in reflected energy followed by a drop as the second cutoff for Steel becomes activated. 
This phenomenon is somewhat easier to observe in figure 3 than in figure 2 because the 
negative phase velocity branch for the AS case occurs in the fourth mode. However, for 
the AS case this negative phase velocity region is present only for a very small range 
of frequencies as opposed to the third S mode where it is far more conspicuous. The 
points /? 2 -« T 2 denote the activation of the third mode and /? 3 ,T 3 the activation of the 
fifth mode. Clearly visible are the windows in the frequency domain where reflection and 
transmission are approximately uniform followed by steep drops in the neighbourhood of 
the cutoffs. It appears from these curves that the transmission of energy is done primarily 
by lowest flexural mode. The higher order AS modes which incorporate most of the 
shearing deformation is poorly transmitted. A similar feature is observed in the AS case 
for other combinations of materials. 

In figures 4 and 5 are shown the symmetric responses for Aluminum to Copper and 
Copper to Aluminum respectively. The properties for Copper are a P-wave speed of 4.17 
Km/sec, an S-wave speed of 2.15 Km/sec and a density of 9.S gm/cc. Copper is a denser 
but softer material than Aluminum. 

An interesting feature of these curves is that the transmission curve for the lowest 
mode is exactly identical whether the transmission is from Copper to Aluminum or vice 
versa. Figures 6 and 7 show that the same is true for the lowest AS mode as well. In 
Appendix B we provide an outline of the proof of this ’reciprocity" for any pair of elastic 
materials. It must be kept in mind that the input excitation consists of the lowest mode 
and hence the ’reciprocity’ is observed only for the lowest transmitted mode. 

Since the cutoffs for Copper and Aluminum are well separated, one can observe the 
phenomena discussed earlier in a somewhat clearer manner. What is particularly remark- 
able in figure 5 is the large reflection (about 60 percent) in the third S mode (dashed 
curve) for Aluminum. In figure 4, however, the transmitted energy is largely in the first 
and second symmetric modes of copper. The blips observed on either the reflection or 
transmission curve are usually stronger at the thickness shear cutoffs as opposed to the 
thickness stretch cutoffs. This is true in transient loading problems as well, where the 
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singularities at the thickness shear cutoffs dominate over their dilatational counterparts. 
In this case, there is a spectral window ranging from about 0.2 to 0.3 Mhz in Copper where 
transmission in the second S mode is quite strong. The blips at f? 4 , T 4 , R 3 , T 3 , i? 5 , T 5 occur 
at the cutoffs of either material. These represent a sudden singe of energy followed by a 
drop or vice versa once again due to the repartitioning of energy. 

In figures 6 and 7 is shown the response to an antisymmetric excitation. As discussed 
earlier there is little transmission in the higher modes. In figure 7 the dotted and dashed 
curves display a strong reflectivity. The same is true for figure 6 where the cutoffs (in 
Aluminum) occur at higher frequencies them those in Copper. 

For the problem of transmission of energy across an infinite boundary, ’mechanical 
impedance’ (the product of density and wave speed) is sometimes used as a measure of 
mismatch between bonded materials. Roughly speaking, this quantity is a qualitative 
measure of one material’s ability to let energy enter or exit into a second material to which 
it is bonded. 

In figures 8 and 9 are shown the results for Aluminum- Glass. These two materials 
are remarkably well matched in terms of their mechanical impedances. For Glass the 
wave speeds chosen are: a P-wave speed of 5.80 Km/sec and an S-wave speed of 3.35 
Km/sec, while the density is taken to be 2.5 gm/cc. Near perfect transmission can be 
observed in figure 8 except near the negative phase velocity region of the third S mode for 
Glass. Almost all the transmitted energy is carried by the fundamental mode for both the 
symmetric and antisymmetric excitations. The reflection is a thousand times smaller than 
transmission with the peak occuring in the negative phase region of the third S mode in 
figure 8 . For the antisymmetric case, the transmission is still better in the fundamental 
mode. The higher modes for both the S and AS cases transmit negligibly tiny fractions of 
energy. 

One of the basic features observed in the foregoing discussion is the interplay between 
the modes in their energy carrying capabilities. As real modes become activated in the 
two materials there are corresponding increments or drops in the energy levels in reflection 
with exactly the opposite being true for transmission. One can view the activation of a 
new mode as an enhancement in the energy transport capabilities of a given plate material. 
A second feature of the results presented is the appearance of windows in the frequency 
domain where a nearly uniform behaviour in energy transmission in certain modes can 
be seen. More work needs to be done to clarify all of the observed features in the plots 
presented here. The goal for future work is to obtain solutions to the transient loading 
problem in bimaterial plates which are of particular relevance to NDE. The technique 
described in this paper can be easily extended to the case of multiple sandwich sections 
which may be of interest in waveguide analysis. Work along these lines is currrently in 
progress. 
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Tables I and II 



Antisyrnm-Refiected Field 


2 Modes 


Mhz 


Mode 1 


Mode 2 


Mode 3 


Mode 4 


0.767 


0.0794 


0.0326 


0.0000 


0.0000 


0.393 


0.1060 


0.0017 


0.0000 


0.0000 


6 Modes 


0.767 


0.3120 


0.0003 


0.0559 


0.1220 


0.393 


0.1510 


0.0058 


0.0000 


0.0000 


10 Modes 


0.767 


0.1900 


0.0181 


0.0369 


0.0438 


0.393 


0.1510 


0.0057 


0.0000 


0.0000 


20 Modes 


0.767 


0.1920 


0.0179 


0.0369 


0.0454 


0.393 


0.1510 


0.0057 


0.0000 
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Appendix A 



For a uniform, isotropic plate there exist both symmetric and antisymmetric Rayleigh- 
Lainb modes which independently satisfy traction free boundary conditions on each surface 
of the plate (y = Hi). The displacement can be written in vector form as follows: 



(1A) u (n) = U (n) (j/)e ,fc " I e- ,t ' t 

where the are either the symmetric or antisymmetric modes. For example, the 

symmetric modes are given by, 



(2A) 



( ik n ch(j n ij) + 6 n B n ch(6 n y)\ 
\Tn^^(TnI/) ik n B n sh(8 n y) J 



The symbols sh and ch refer to the hyberbolic sine and cosine respectively. The quantity k n 
which is related to the horizontal component of the wavenumber satisfies a transcendental 
equation called the Rayleigh-Lamb frequency equation which for symmetric modes can be 
written, 



(3A) (2 k’l - K 2 ) 2 ch(^ n )sh(6 n ) - 4k 2 n ^ n 6 n sh(-y n )ch(8 n ) = 0 

where 



(4A) 

(5A) 



6 n =Hi- !<■)'/■ 



(CA) B„ = (2 k;, - K-)ch(-,„)/2ik„6„cli(6„). 

In the above, k and K are the compression and shear wave numbers respectively. The 
corresponding stress vector for the n th eigenmode is: 



( 7 A ) 






where again the can be either the stress resulting from the symmetric or antisymmetric 
displacements. For the symmetric case it is given by: 



(SA) 



c(n)/, \ _ ( i.~Kn + I\ 2 )ch(~, n y) + 2 ik n B„8 n ch(6 n y)\ 
2 ik nln sh( ln y) + (2 kl - K*)B n sh(6 n y) J 



The bi-orthogonality condition for the in-plane problem can now be simply stated as 



(9A) ,, - Vi""S%>)dy = ( 0 * " 

J-h • [ Jm m = n 

Here the subscripted forms of S and U represent components of each vector. Notice that 
the biorthogonality relation involves an integration over the thickness of the plate. 
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Appendix B 



Following is a short argument proving the reciprocity between transmitted fundamen- 
tal modes of a bimaterial plate noted in the results section of the paper. 

In the absence of body forces, the reciprocity theorem (Achenbach et al (1982)) 
between two elastic states (A and B) can be written: 

(Bl) J (ufr,j - ufr? )rijdA = 0. 

For the 2D plane strain bimaterial plate problem with plate thickness 2h, the above integral 
becomes: 



(B2) 




V 4 

i l ij 




r B 

O 






where the interface between the two plates is at x = 0. 

For elastic states A and B consider transmitted and reflected fields (expressed in terms 
of Ravleigh-Lamb eigenmodes of the appropriate plate) resulting from incident displace- 
ments 



(B3) ui = 

(B4) = 

where J] and Q\ are given in equations (12a-b). By substitution of the above incident and 
resulting reflected and transmitted fields into equation (B2) one finds that the coefficients 
of the leading terms in the transmitted fields of each elastic state are equal. 

It should be noted that if the incident field involved additional modes the reciprocity 
relationship would involve a combination of the corresponding transmitted modes. 
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